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Abstract 

Plastino and Curado [Phys. Rev. E 72, 047103 (2005)] recently determined the equi- 
librium probability distribution for the canonical ensemble using only phenomeno- 
logical thermodynamical laws as an alternative to the entropy maximization pro- 
cedure of Jaynes. In the current paper we present another alternative derivation of 
the canonical equilibrium probability distribution, which is based on the definition 
of the Helmholtz free energy (and its being constant at the equilibrium) and the as- 
sumption of the uniqueness of the equilibrium probability distribution. Noting that 
this particular derivation is applicable for all trace- form entropies, we also apply it 
to the Tsallis entropy showing that the Tsallis entropy yields genuine inverse power 
laws. 
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1 Introduction 

The statistical mechanics tries to explain the thermodynamic behavior of 
macroscopic systems by relying on the microscopic constituents of the physical 
system. In this context, it is of utmost importance to obtain the equilibrium 
distribution associated with a particular entropy structure corresponding to 
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the thermodynamic properties of the system. A well-known approach to obtain 
the equilibrium distribution in statistical mechanics is the entropy maximiza- 
tion procedure of Jaynes (JF) [T]. This maximization procedure in turn stems 
from the second law of thermodynamics, since entropy attains its maximum 
value only at equilibrium. 

On the other hand, it is necessary to be able to obtain the equilibrium distribu- 
tions purely on thermodynamic grounds by also including some microscopic 
ingredients such as microstates in order to have a comparison with the en- 
tropy maximization of Jaynes. This necessity is twofold: First of all, it is not 
straightforward to conclude the validity of second law for generalized entropies 
[2|5] . Second, it is apparent that even ordinary statistical second law requires 
a deep and detailed study [1]. 

The generalized entropies do not change the inherent dynamics so that the first 
law of thermodynamics (FLT) and its solutions must be certainly preserved 
i.e., they must be form invariant independent of the entropy structure one 
adopts [3]. This is also required, since the first law of thermodynamics is the 
conservation of energy [B]. The second law of thermodynamics too must be 
preserved in these generalized schemes in order to construct a link between 
information-theoretic entropy measures and the statistical mechanics. In short, 
then, the generalized entropy measures must conform to the thermodynamical 
laws in order to be considered useful in statistical mechanics, instead of being 
solely useful for information theoretical approaches. 

The paper is organized as follows. The next section outlines how to obtain 
canonical equilibrium probability distributions through the requirement of a 
constant Helmholtz potential F and the assumption of the uniqueness of the 
aforementioned distribution. This approach is then extended to all trace-form 
entropies and applied to the Tsallis entropy measure [7]. Conclusions are pre- 
sented in the last section. 



2 Thermodynamics and Equilibrium Distributions 

The first law of thermodynamics for reversible processes reads 



where U is the internal energy, S is the entropy and T is the temperature, 
while P and V denote pressure and volume, respectively. The term /i is the 
chemical potential associated with the change in the number of particles A^. 

A particular case of great interest in thermodynamics is when U = U{S, V, N) 
is a homogeneous function of first degree [B]. In this case, Eq. ([T]) can be 



dU = TdS -PdV +ijdN, 
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explicitly solved, independent from the analytical expressions of T, P and /i 
[8], yielding 

U = TS - PV + fiN . (2) 

The equation above is referred to as the Euler equation of thermodynamics. In 
terms of the free energy F = F{S, V, N) (also called the Helmholtz potential), 
the above equation can simply be rewritten as 

F = U-TS. (3) 

Keeping in mind that F = —PV + fiN, Eq. ([3]) is equivalent to Eq. ([2]). This 
definition of Helmholtz potential F is also in accordance with the rest of the 
thermodynamical potentials. For example, the Gibbs potential G is equal to 
F + PV so that one obtains the well known relation G = fiN. We can further 
rewrite Eq. ([3]) as 

S = I3U-(3F. (4) 

where 1/(3 = T. 

Up to this point, all these equations contain nothing related to the statistical 
nature of the physical system under investigation. In order to introduce in 
Eq. (j3]) statistics, one needs some statistical ingredients. We then consider 
the very definition of Boltzmann- Gibbs (BG) statistical entropy, which reads 
(setting /cb = 1) 

Sbg = J2Pi^^(^/Pi) ^ (5) 

i 

and the definition of the internal energy in terms of the microscopic proba- 
bilities i.e., U = J^iPi^i, where pi is the probability of the corresponding ith 
microstate and Ei is the microstate energies. 

Then, when the Helmholtz potential is constant, Eq. (j4]) can be rewritten as 
^pan(lM) =(3J2p^e, -f3j2p.F, (6) 

i i i 

where we also assumed probability normalization as usual. Eq. ([6]) naturally 
yields to the relation below 



In (l/pi) - /3ei + f3F 



0. (7) 



However, it is not trivial to obtain the equilibrium probability distribution 
p* from the equation above, since whole summation might be equal to zero, 
and not only one singled out term within the summation. On the other hand, 
assuming the existence of the thermodynamic state of the stationary equilib- 
rium, we can invoke the uniqueness of the aforementioned state and therefore 
of the corresponding equilibrium distribution. Due to the former assumption 
i.e., the existence of the equilibrium state, we have p* ^ 0, since the condition 
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p* = contradicts with the existence of an equihbrium state. The unique- 
ness of the thermodynamic state of equihbrium and hence of the equihbrium 
distribution imphes the trivial solution to Eq. ([7]), equating the term in the 
brackets to zero i.e., 

In - + /3F = 0. (8) 

It is then obvious that the equation above yields an equilibrium distribution 
and this is all one needs, since the equilibrium distribution (and the state of 
equilibrium) is unique, so that having one solution of Eq. ([7]) is tantamount 
to obtaining the only solution indeed. Therefore, we finally have the following 
equilibrium distribution 

p^ = exp{-P{e,-F)), (9) 

which is the canonical equilibrium distribution. Moreover, note that identi- 
fying e~^^ as the canonical partition function Z, one obtains another well 
known relation i.e., F = —TlnZ (remember that we set Boltzmann constant 
/cb to unity). We observe that the canonical distribution p* satisfies Eq. (jlj) 
and the FIT in Eq. (II]). 

Some remarks are in order: although mathematically simple, the above calcu- 
lations are novel in the sense that it does not make use of ordinary entropy 
maximization procedure put forth by Jaynes. Apart from the definitions of 
entropy and internal energy in terms of microscopic constituents, it relies 
solely on the definition of Helmholtz free energy and the existence and the 
uniqueness of the thermodynamic equilibrium state. Moreover, obtained in 
this manner, one sees that the canonical BG-equilibrium distribution, apart 
from the multiplicative factor /3, is of exponential form, where the exponent is 
formed by shifting the microstate energies Si by the constant (and minimum 
at the canonical equilibrium state) Helmholtz free energy F. 

The above calculation can easily be extended to generalized entropy measures, 
since these definitions are information-theoretic and not supposed to violate 
any of the ingredients used in the above derivation. In order to extend the 
formalism above to generalized trace-form entropies, we study the following 
statistical representation of the thermodynamic entropy i.e., S = (A(l/pj))j.j^, 
where ((■ ■ ■ )j)ijj^ = J2iPi{- ■ ■ )« is the arithmetic expectation value of the quan- 
tity (■ ■ ■ )i, Pi are the microstate probabilities, and assuming that the internal 
energy is given statistically as f/ = {Bi)i^^, where et are the microstate energies. 
Similar to Eq. we now write 

(A(1M))h„ = mur. - (/^^)iin =^ - Pe, + = . 

(10) 

Note that A{l/pi) is a generic function, which can be any monotonic increasing 
function compatible with the linear averaging scheme, be it ordinary or gen- 
eralized logarithm. ^From the relation above, we can obtain the equilibrium 
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probability distribution pi, which is compatible with FLT, namely 



A-i(/3(£,-F)) 



Z 



S = P{U-F) = Y.PWe^ ©A KZ)] , (lib) 

i 

where x ©a V ■= A(A~^(a;)A~-'^(?/)) {x,y > 0) and ©A^in + is the ordinary 
addition operation. Z in Eq. f lllap represents the partition function, defined 
implicitly through Eq. f lllb[) . As can be seen, its analytical expression depends 
on the function set {A, A~^}. If we now substitute the probability distribution 
given by Eq. flllap into the left hand side of Eq. (jlj), it reproduces indeed the 
Euler equation. 



2.1 Tsallis entropy 



An entropy measure of current interest is the Tsallis entropy [711^ . This entropy 
can be constructed by choosing 

A(g, x) = \n,{x) = ^JZ^ ' ^) = e^P9(^) = [l + (1 " ^ (12) 

where = max{0,x}, so that the structure S tends to the following q- 
entropy 

^. = Ep^K(iM) = ^f!~^ (13) 

The equilibrium probability distribution, according to Eqs. f lllap and flllbP . 
takes the form 

— = expg {Pisi-F)) = Zg expgi^Ei) , with = exp^ 

(14) 

We note that this equilibrium distribution structure is exactly the one obtained 
through the entropy maximization procedure of the Tsallis entropy with linear 
constraints [TfTO] . It is also worth noting that this method apparently does 
not apply to non-trace form entropies. 



3 Conclusions 



In this work, we derived the well-known canonical equilibrium probability dis- 
tribution by using the microscopic definitions of entropy and internal energy 
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together with the definition of the Helmholtz free energy (and its being con- 
stant at the equihbrium) and the existence/uniqueness of the state of equihb- 
rium. Providing a generahzation of this derivation to all trace-form entropies, 
we also obtained the stationary equilibrium probability distribution of the 
Tsallis entropy which is found to be an inverse power law distribution. This 
work can be considered as an alternative to both the entropy maximization 
procedure of Jaynes [1] and recently proposed procedure of Plastino et al. 
[TT|12fl3|14|15j based on the phenomenological laws of thermodynamics. 
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